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ABSTRACT 

The HEAO-1 experiment A-2 has performed a complete X-ray survey of the 
8.2 steradlans of the sky at |b| > 20® down to a limiting sensitivity of < 3.1 
X 10”^^ ergs/cm^ sec In the 2-10 keV band. Of the 85 detected sources 
(excluding the LMG and SMC sources) 17 have been Identified with galactic 
objects, 61 have bean identified with extragalactic objects and 7 remain 
unidentified. The log N - log S relation for the non-galactlc objects is well 
fit by the Euclidean relationship. We have used the X-ray spectra of these 
objects to construct log N - log S in physical units. The complete sample of 
identified sources has been used to construct X-ray luminosity functions, 
using the absolute maximum likelihood method, for clusters of galaxies and 
active galactic nuclei. 

Keywords: X-Ray Sources, Luminosity Function, Cosmic X-Ray Background 
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I. INTRODUCTION 

Th« HEAO-1 aateliita experiment A-2 (Bothichlld et al. 1979) with Its 
extended energy range, complete sky coverage, low and stable internal 
background and moderate spatial resolution has enabled us to create a complete 
catalog of X-ray sources at galactic latitudes |b| > 20° down to a limiting 
fsnsitivity of 3.1 x 10"^^ ergs/cm^ sec in the 2-10 keV band. Recent 
identifications of these sources by modulation collimator experiments on 
HEAO-1 and SAS-3 as well as imaging detectors on HEAO-2 has resulted in 
certain identlf tcations of all sources of flux ^ 4.0 x 10"^^ ergs/cm^ sec, 
pending confirmation of two clusters and NGC 7172, and reasonable 
Identifications for 78 out of the 85 (92%) sources in the sample. All but 9 
of these identifications are extremely likely or certain. This identification 
ratio for the extragalactlc sources compares to identification of 45 out of 67 
(67%) sources in the sample of Warwick and pye (1979). 

The completeness of this sample enables construction of the 
number-intensity distribution (log N - log S) for X-ray sources as well as 
developing X-ray luminosity functions for clusters of galaxies and active 
galactic nuclei. In addition the body of X-ray spectral data returned by A-2 
allows us to cast the log N - log S distribution in absolute rather then 
instrument dependent units which enables comparison with the log N - log S 
relation in different X-ray energy bands (cf. Giacconi et al. 1979). 

Analysis of this data shows that the source counts are well fit by a 
"Euclidean" law with 

dN -5/2 , -1 -1 

^ - 16.5 S (R15 cts/sec) sr 

consistent with previous results despite the quite different samples (Warwick 
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and Py« 1978; schwartfe 1979)* Th« luaihoalty functions at* wall fit by powar 
lau raprasantations with 

dN < , _ ,„-7 -2,15 ,,„44 . ,-1 „ -3 

— ■ 3,5 10 (10 arg/sao) Npc 

44 

for clustats of gaiaxies* and 

^ ^ ,-2.75 ,,^44 , ,-1 „ -3 

•3T ■ 2.7 X 10 44 (10 erg/aac) Mpc 

44 

for active galactic nucleii almilar to previous results (McKee et al. 1980} 
Pye and Warwick 1979). Integration of the luminosity functions over 
the < 10^^*® - 10^® erg/sec range within which they are well determined 
results in estimates of the contribution of clusters and active galactic 
nuclei to the Integral 2-10 keV unresolved X-ray background of 5 4% for 
clusters and < 20% for active galaxies. Using these luminosity functions, 
with no evolution, wo estimate that ^ 30% of the sources seen in the Einstein 
Observatory deep survey (Giacconi et al. 1979) should be relatively low 
luminosity (L < 1 x 10^^ erg/sec) nearby (z ^ 0.5) objects. 

II. DATA ANALYSIS AND SOURCE SELECTION 
The HEAO-1 A-2 experiment, described in detail by Rothschild et al. 
(1979), provided two independent, low background, high sensitivity surveys of 
the entire sky six months apart. We have analyzed the A-2 data in order to 
obtain a complete flux limited sample of extragalactic X-ray sources. The 
region between -20® and +20® in galactic latitude has been excluded to 
minimize contamination from galactic sources. A circle of 6 degrees radius 
around the LMC sources has been also excluded to prevent confusion problems. 
Therefore, we remain with 65.5% of the sky (8.23 Star). The statistical 
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■ignific4nc« of tho «xiit«nc« of tha •ouroaa .la taatad Ijy datarminln^ tha 
dacraase in x whan tha na«r adurca la addad to tha modalt All aouroaa in tha 
aampla giva a dacraaaa in x «t laaat 30 • Tha probability of having, by 
chanca, a dacraaaa of 30 in x wibh two dagraaa of fraadom (acan angla and 
intensity) la 3 x 10~^t Thia probability la almost tha aama aa the one 
aaaooiatad with a deviation of So in a Gauaaian diatribution (6 x lO”^). 
iherefore, we can alao state that the lowest statistical significance for the 
existence of the sources included in our sample is So, as required by the 
maximum likelihood methods we use to determine the log N - log S parameters 
(see Section IV- 1). Taking into account this statistical significance 
requirement we estimated the completeness level of the first and the second 
scan as 1.25 and 1.8 RIS counts/sec respectively, see Figure 1. One R15 
count/sec ^ 2.17 erg cm”^ sec“^ in the 2-10 keV energy band for a power law 
spectrum with photon index 1.65. R15 is a counting rate der.ived using the 

1.5° X 3° FWHM fields of view of the second layer of the argon counter and 
both layers of a xenon counter. This combination has a FWHM for the quantum 
efficiency from < 3 to < 17 keV (Marshall et al. 1979). 

The second pass is less sensitive on average, because much more time was 
spent in pointing at sources. We shall be more concerned with the first pass 
data in deriving best fit parameters and use the second pass ones mostly as an 
independent confirmation. 

III. OBSERVATIONS 

A. The Sample 

Table 1 contains all the relevant data for the 68 sources either 
brighter than 1.25 R15 c/s in the first scan wh.ich corresponds to days 248-437 
of 1977, or brighter than 1.8 in the second scan, days 73-254 of 1978. Source 
names are listed in column 1. Column 2 contains previous catalog names. 

I 
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Firat pass £lux(»s and 1» arrors ara In column 3# Mhila tha saoond pass onas 
ara in col'Omn 4 . Soma fluxas may dlffar slightly from pravlously raportad 
results, as different procaduras have bean usadi a.g., in tha recant paper by 
McKee at al. (1980) fluxes have been obtained fixing tha X-ray position at the 
optical position, instead hare we have used the beat fit X-ray position to 
derive the flux. Available identifications are listed in column 5. The type 
of object is in column 6. One • in column 7 indicates firm identifications 
(i.e. as provided by tha SAS-3 ot HEAO-1 modulation collimator/ or the 
Einstein X-ray telescope), two * indicates possible identification consistent 
with larger error boxes. Redshift values and references are given in column 
8. Spectral information is now available for more than half of our sources 
(Mushofczky et al. 1980/ Worrall et al. 1980> Hushotzky 1979i Holt'l980» Boldt 
1980), we quote in column 9 conversion factors between R15 counts/s and ergs 
s“^. Wi.en spectral information is lacking we assumed a 6 keV thermal 
bremsstrahlung spectrum for all sources identified with clusters and a 1.65 
photon index power law for all sources identified with active galaxies. An 
average conversion factor value of 2.5 x 10“^^ ergs cm”^ b"Vr 15 counts s"^ 
was assumed for the few unidentified sources. Columns 10 and 11 contain the 
first and second pass luminosities in units of 10^^ erg s“^ calculated for 
K 50 kra/s/Mpc and = 0.5. Column 12 contains notes. 

B. Classes of Sources 

Sixty of the 82 sources brighter than 1.25 counts s"^ in the first scan 
and not definitely associated with galactic objects have been associated with 
extragalactic objects. Only 7 remain unidentified at the present time. These 
60 identified sources subdivide almost equally between clusters of galaxies 
(30) and single galaxies (30). Most of the 30 galaxies are Seyfert galaxies 
of class 1 or 2, but we have also 1 QSO (3C 273), 4 BL Lac objects, and 1 
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"nomal** galaicy (NGC 71,72) • Not* that N31 and tha Magallanic Cloud aouroaa 
ata not Inoludad in our uxtragaiatlc Bampla bacauaa thay rapreaant a local 
inhomoganalty aa part oC tha local group oC galaxlaa. Tabla 2 llata tha 17 
high galactic latituda aourcaa not included in our axtragalactic aampla 
bacauaa thay have been identified with galactic or "local** objecte« The 
second pass sample contains only 37 sources brighter than 1.8 Rl5 counts/sec, 
all but one identified. The source classification is consistent with the 
first pass. Assuming Poisson errors, clusvers contribute 50 ± 9% of the 
identified sources in the first pass and 61 ± 13% in the second. Galaxies 
contribute 50 ± 9% in the first scan and 39 ± 10% in the second. 

C. New Sources and Sample Completeness 

H0328t025 and H0917-075 are the only entirely new sources in Table 1. 
Figure 2 shows their error boxes. All the other sources in Table 1 have been 
listed somewhere else before. The Improvement in our sample, as compared to 
previously reported ones, is due to a better rejection of non-extragalactic 
sources, made possible by the recent identifications, and to a uniform sky 
coverage to a relatively low limiting flux. 

As the instrument has a fairly large (1.5® x 3.0®) angular resolution 

the possibility of source confusion must be considered. The total area of the 

sky Included in this survey is approximately 2.7 x 10^ square degrees, 

therefore there are about 6 x lO^ independent positions on it. As the high 

galactic latitude X-ray sources bright enough to give confusion problems at 

our sensitivity level cannot be more than a hundred using the log N-log S 

relation derived later (taking into account also the possibility that two 

weaker sources can simulate a source bright enough to be included in our list) 

we therefore expect negligible confusion. That is using ~ « 16.5 s”^*® there 

qs 

are roughly 65 resolution elements per source, of S > 1.25 cts, well above the 
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confusion l«v«l of 25 boan aroas p«r aourco oftan quotad in tha litaraturat 
In addition tha uniform aky coyaraga at tha choaan aanaitivity lavala providad 
hi' this axparimant and tha availability of two indapandant aata of data for 
croaa-chacking purpoaaa support our confidanca in tha conplatanaas of our 
sampla* 

p. Space Distribution of Sourcas 

since the pioneering work of DeVaucoulaurs (1958) much attention has 
been devoted to finding evidence of a superclustar centered in the Virgo 
cluster of galaxiest We plotted the positions of our sources in supergalactlc 
coordinates looking for some kind of anisotropy • Figure 3 shows the 1st pass 
sample. Obviously > no anisotropy is observed as most sources lie beyond the 
supercluster. If we restrict our attention to the 12 sources with redshifts 
less than .01 (in boxes in Figure 3), we see that 9 are in tha center region 
of the supercluster while 3 are in the anticenter and that all but one have 
supergalactlc latitude less than 30 degrees in absolute value This result, 
which is significant at the few percent level, suggests that close X-ray 
galaxies may lie preferentially in the supergalactlc plane. But no conclusion 
can safely be Made from such a small number of objects at present. 

IV. THE NUMBER-FLUX FUNCTION 
The usual power law form 

“•(X • i • 1 

N(S) ■ KS (R15 counts/sec) 'or (1) 

has been assumed for the nUmber-flUx relation. The various methods applied to 
estimate the coefficient K and the differential exponent a as well as to 
evaluate the goodness of the fit are outlined in the next section. 

A. Statistical Methods 
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1. NaxiJRURt Llkalihood 

Crawford, Jauncay ai..d Murdoch applied the maxiiMiiK likelihood Mthod 
to unbinned data in order to eatimate the elope off the nuMber«'lflux relation of 
radio aourcea* In the firat paper (Crawford et al. (1970) a aolution ia 
worked out for error free data, in the aecond paper (Murdoch et al. 1973) the 
method is extended to Include eirrora on the measured fluxes. Numerical 
corrections to the error free answers v^are calculated for the special case of 
Gaussian distributed errors. In the aame paper it was pointed out that a 
minimum signal- to-noise ratio o| five is required so that the uncertainty in 
the correction factor due to weaker sources does not dominate the correction 
itself. This is why we excluded from our sample sources with statistical 
significance less than So. Xn both tmpers i and XI the Kolmogorov-Smirnov 
teat (here after: K-S test) was suggested to evaluate the goodness of the fit 
obtained. Xn the remainder of this paper we will refer to this method as to 
the Maximum Likelihood (ML) method. 

2. Absolute Maximum Likelihoed 

The ML method assumes the same underlying error distribution for 
all the souces in the sample, i.e. it assigns the same lo error to all the 
sources. As we deal with sources of greater than So significance the error 
assumed is one fifth of the minimum flux in the sample, or .25 R15 counts s"^ 
in the first scan and .36 R15 counts s"^ in the second. Table 1 shows that 
these values are not very far from the actual errors. However Lightman et al. 
(1980) have developed a refinement of the ML method in connection with the K-S 
test capable of handling sources with their own experimental error. Following 
those authors we will call this new statistical method the "Absolute Maximum 
Likelihood" (AML) method. Lightman et al. (1980) worked out the AML method on 
general grounds and then applied it to the evaluation of globular cluster 
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X-ray aourqa masaaa< h» thia is ths (first application of tha AML mathod to 
tha numbar flux funotioni wa giya a short outline of tha isathod balow. 

Assuming tha form (1) for tha numbar-flux ralatlpn and a @aussian 

form for tha arror distrlbjition of tha msasurad fluxas wa avaluatad 

numarically tha integral probabilitias (a) as 


/ *-dF X" d$ N(S,o) p(P,o.,S) 
§ . -Jaifl 52 

^ ® X* ds N(s,0) p(F,o. ,s) 

*^min CO 


( 2 ) 


where S is tha true flux, and are the measure*', flux and error of the 
i-th source. is a cutoff value used to avoid the apparent divergence at F 

■ 0. As in Murdoch et al. Cl‘>73) the particular choice of the cutoff value 
does not affect the value of the integral as long as the statistical 
significance of the sources is at least 5c. ^mln sensitivity limit 

of the sample. For every assumed a we computed the Fj^(o) for all the 
sources. The Pj^(a) should be uniformly distributed between 0 and 1. 

Following Lightman et al. we evaluated the maximum deviation from the uniform 
distribution: 


D (a) 
max 


» max [D.(o)] 
i«l,N 


- max ( |Pj^(o) " N * ^ 

i"sl,N 


(3) 


where N is the number of sources in the sample and the Pj^(a) have been sorted 

in ascending order. Then we calculate the probability (a)) of 

observing deviations greater than from the formula for the K-S 

statistic given by Birnbaum and Tingey (1951)n The (a,P(D (a)) function is 

max 

then plotted. The best fit value of a is the one corresponding to the maximum 
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v«Lu« of th« distribution. Obviously iwjst bs grsstsr 

than soms mininuin vslus (ssy 10%) In ordsr to socspt th# Modsl. Ths rsngs 
in a fiiif rssults givsn bslow on a am svslustsd Croii ths vsluss 
«nd Ojj of a, which raduos P(Cj,|,^(o)) to Pk/^){'^2. 

3.Chi-Squars 

Both ths ML and ths AML msthods ars indspsndsnt of ths 
cosfficisnt ic of ths numbsr-flux rslation, as K Is lost In normalizing ths 
probabilitisa. Thsrsfors, wa ussd ths mathod to dstsrmins ic. Bins with 
squal sxpsctsd numbsr of sourcss for a ■ 2.5 havs bssn ussd for 
ths calculations. Of course# in calculating confidsncs bounds# we havs 
assumed that the functional form of the distribution is the "true** one* If 
better data later shows that this Is not true our confidence values are not 
applicable. 

V. LOG N - LOG S RESULTS 

The ML method applied to the 60 non-galactlc sources brighter than 1.25 
R15 counts s’^ in the first pass gives (in this section wo use R15 counts s"' 
as the unit) 

0 - 2.67 ± .23 

with a goodness of fit probability (evaluated using the KS test) of 39.5 
percent. 

For the 37 non-galactlc sources brighter than 1.8 R15 counts s”i in the 
second pass the ML result is 


a « 2.74 i .32 
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nith A probability of 17*5 poroont. 

AML cooultm aro 

O " 2>63 1 <2 

in tha first paaa« aao Pigura 4a# and 
0 - 2.74 ± .22 

in tha second# sea Plgura 4h. 

*rhe 68 and 95 parcant probability contours for the 1st pass values of k 

2 2 
and a evaluated with tha x method are plotted in Figure 4c Tha x baa^ ^ib 

values and 1o errors for the number-flux function parameters are 

a « 2.72“^^Jg (4) 

m> 2#6 fli*! 

K " 20_^^|q (R15 counts/sac) er . 

The differential number-flux data as well as the bast fit function 

N(S) “ 20 (R15 counts/sec)”^ sr"^ 

2 

are plotted In Figure 5; the y value of the fir Is 2.79 for 6 degrees of 
freedom# corresponding to a probability F( > X ) ^ limited size of 

the second scan sample does not allow a good estimate of the probability but 
the results are consistent with the first pass ones. 

C. Number Flux Relation in Physical Units 

Using the conversions factors listed in column (9) of Table 1 we can 
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•xpreas th« fluxes In ergs cm"^ s"^ end evaluate the number- flux relation 
accordingly* Conversion factors range approximately from 2*0 x 10"^^ to 2*9 x 
10"^^ ergs cm’^ ,“1 (rib counts s"'')'^' highest values referring to soft 
spectra sources whose emission peak lies below our Instrument energy window. 

As a consequence of the different conversions factors the completeness level 
of the samples when fluxes are Ln ergs cm~^ s~'^ is equal to the former 
completeness level in R15 counts s"^ times the maximum conversion factor: 
that is < 3.6 X 10"^^ ergs cm"^ s“^ for the first pass and 5.2 x 10"^^ 
ergs/cm^ sec in the second pass. The 1st scan sample with this flux 
restriction contains 51 sources: 25 clusters# 22 "galaxies'* and 4 

unidentified sources. The best fit values and 1 a errors for the number-flux 
function parameters obtained with the three methods agree with 


a ± 2.85 ± .3 

7 /r / -2 -IvO-l -1 

K < (5.65,- _) X 10 (ergs cm s ) sr 

"" ▼ 1 s 


(5) 


The 32 second scan sources brighter than 5.2 x 10"^^ ergs cmT^ s~^ give us a 
best fit of slightly steeper slope a 3.1 ± .4. 

VI. DISCUSSION OF THE RESULTS 

All the first pass samples# whether fluxes are e:q>ressed in R15 counts 
s"^ or in ergs cm“^ s"^ are consistent with the five halves Euclidean slope 
(see Figures 6 and 7). The slight preference for a steeper than Eucledian 
slope is due to the distribution of the brightest few sources in calculating 
the likelihood functions. It is these sources that are most sensitive to 
changes in a by virtue of the relatlveljf small statistical error in their 
measured intensity. Our Euclidean best fit for the 1st pass data is 


N(s) - counti/«)“%r'’' 

with « d*gic««» oC frandomi p(x^ > 5.5) < 60%. Th« ML 

probabiXity for a- 2.5 is 42.4%. Asaumln? an av«ra 9 « convarsion factor of 2.4 
X 10“^^ argt cm"28‘“^ (R15 count! a)“^ tha ralatlon (4) bacoinea 

N(S) < <1*®^1?25^ ** cm"^a“^"^ ar"^ 

in agraament with tha exact raauXt 

N(S) < (2.2^*'^) X Xo“^^ s"^’^ (arga cm"^a“S~^ ar^^ 


obtained from the lat paaa complete sample for fluxes in ergs cm“^s“’’ and 
using the conversion factors in Table 1. 

VII. COMPARISON WITH PREVIOUS EXPERIMENTS 
Both the Uhuru (Schwartz 1979) and Ariel 5 data (Warwic)c and Pye 1978) 
gave a flux-numbar function consistent with the Euclidean model. Their best 
fit values for the coefficient < with o ■ 2.5 and S in Rl5 counts a"^ are 
respectively 

(c « 16.5 ± 3.9 uaing 1 Uhuru ct/s »• 1.0 R15 ct/sec 

and 

X "• 15.8 ± 4.2 using 1 Ariel-5 ct/sec •» 2.12 R15 ct/sec 

in agreement with our results at the 1o level. These conversion factors 
assume a mean R15 conversion factor of 2.4 x 10"'^’ ergs/sec, 1 UhvJru ct/sec « 
2.4 X lO""*^ erg/sec, and 1 Ariel-5 count/sec « 5,1 x 10"'*^ erg/seo. If we use 
the calibration of Marshall et al. (1979) appropriate for the active galaxies 

of 1 R15 ct/sec « 2.17 x 10"^^ erg/cm^ sec, we find t<ohuru ^ 
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tJ-opt of WACwiok and py« of 2.7 i .2 is aXso 

consiatsnt with out tsiuit. 

vm. LUMINOSITV PUNCKON 

h* Msthod and Data Baas 

Many authota (Sohwartie 1970| Moiiatdy I978i MoKaa at al. i980| Blvia at 
al. i977> Pya and Watwick 1979» Tananbauw at ai. i978» Boidt 1980) hava 
tacantiy considacad t)m ptohlam of avaluating the X-ray luminosity functions 
for different olassas of sources principally^ olustars of galaxies and active 
galaxies. Ml of them with the exception of Pya and Warwick l»ad to raly upon 
optical data to select complete samples, wa present hare X-ray luminosity 
functions evaluated from X-ray flux ll?p,ited samples, as we remain with a few 
unidentified sources i our results have soma uncertainty* but wa believe that 
the residual incompleteness should not be very important. 

1 . the Samples 

The first pass sample of clusters of galaxies contains 3Q ob;)eQts. The 
second pass ona includes 22 sources. Thirty •’galaxies*’ are observed in the 
first pass* but wa excluv*,e from our sample the QSO 3C273* the 4 dl Lac 
objects* tile peculiar galaxy M82 and the ’’normal** galaxy NGC 7172 as they are 
not liomogeneous with t)»e bulk of the sample whicli consists of Seyfart 
galaxies. T)>erafora wo remain witl\ 23 active galactic nuclei. The second 
pass sample contains only 12 objects (after excluding 30 273 and PKS 
2155-304). 

The completeness of the sample is checked using the Schmidt > test 

and wl.th a X-S test on the distribution of the as suggested by Xvni and 

Bahcall (1980). The resulhs are listed in Table 3. 


TABLE 3 
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* OBJECTS 


K<'S TEST 

CLASS OF OBJECTS 

SCAN » 

IN SAMPLE 

> 

V 

POd) % 

Clusters of Galaxies 

1 

30 

.4711.054 

18.1 

Clusters of Gslaxies 

2 

22 

.5521.062 

11.8 

Active Galaxies 

1 

24 

.5231.059 

50.4 

Active Galaxies 

2 

12 

.5571.083 

56.8 


The 10 error quoted for <V/Vj^> ia the fomuiX error 1//12 n 7 where N is the 
number of objects In the sample (aee Avni and Bahoall)* All the 4 samples 
meet the requirements of the tests* However^ we expect a sm^^';!.! degree of 
Incompleteness due to the unidentified sources. 

2. Methods of Analysis 

Of the three methods outlines in Sec (IV-A) only the AML is suited for 

the determination of the luminosity function parameters. The relatively small 

2 

sizes of the samples do not allow an efficient use of the x square method or 
of any other binned method. Moreover the ML method in the form developed by 
Crawford, Jauncy and Murdoch cannot be used because of its assumptions of a 
single underlying error distribution. This last hypothesis was reasonably 
satisfied by the flux data in the evaluation of the log N log S parameters, as 
we already pointed out, but is not satisfied at all by the luminosity data, as 
the errors are proportional to the square of the redshift of the sources: 



( 6 ) 


On the contrary the AML method is well suited for the task. The description 
of Section IV-A still applies. However, instead of calculating the 
probabilities of eq (2) we evaluated the probabilities: 
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— (7) 




Eq> (7) gives the integral normalized probabilities of observing a source with 

measured luminosity less or equal to assuming a Gaussian etror 

distribution with standard deviation 0 , and for the differential luminosity 

^i 

function the form f(lj,q) where L is the true luminosity and q represents che 
functional parameters to be determined. and are the lower and upper 

boundaries of the luminosity function. is the maximum volume at which 

one could detect the source and depends on the sensitivity limit of the 
sample. For a source of luminosity L in a sample of minimum sensitivity F^^^^ 
the maximum visibility volume is proportional to 

Note that Eq. (7) does not take in account errors on the redshift 
z. 1*he AML method can determine the form of the luminosity function but not 
its absolute value. Therefore we have used a least squares fit to the 
unbinned data to evaluate the multiplicative coefficient. 

B. Results 

1. Clusters of Galaxies Luminosity Function 

We considered t«K> different forms for the luminosity function: the 

power law form 


f(L) - KL 


and the exponentfial form 
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£( 1 *) « 


between the minimum (l*44n,j_n) and maximum (1*44^^) observed luminosities, 
expressed in units of 10^^ ergs sec“^. The normalization for a power law 
luminosity function scales as Hq~^> 

Clusters of Galaxies 

Figure 9 represents the AML probabilities for the slope of the cluster 
of galaxies power law luminosity function. The 1st pass best fit values for 
the power law parameters are 

Y = 2.15_ 

K» (3.5 ± 1.1) X 10 ^ (10^^ erg/s) ^ ^ Mpc 


K has been evaluated with the least squares method. The error on K has been 
determined by letting y assume the 1o extreme values of 2. 03 and 2.32. Figure 
8a gives a binned representation of the data with the best fit luminosity 
function. Each bin contains three sources, except for the highest luminosity 
bin which contains five. The second pass results are 


Y 


2.13 


+ .16 

-.24 


-7 44 Y-1 -3 

K = (3.8 i 2) X 10 (10 ergs/s)' Mpc 

Figure 8b give the binned representation. The minimum luminosity object in 

both the 1st and the 2nd pass at 2.4 x 10^^ (ergs/s) is the Virgo cluster. 

The highest luminosity cluster is Abell 2142 winh 2.8 x 10^® (ergs/s). 

The exponential form of the luminosity function has also been 

considered, but the quality of the fit is poorer, see Figure 10. 


18 


M th« Vir 9 o CIUBtar of has a "local" chacaotar, w« avaluatad 

tha cluatar of ^laxiaa luninoaity function without tha Virgo cluatar. Tha 
1st paaa satnpla la raducad to 29 aourcaa, tha maan ^ 0>05S and 

cha K-S taxt on tha uniformity of tha dlatrlbutlon glvaa a probability 

of 24. 7% • The 2nd paaa aampla contains 21 aourcaa, tha maan >la 0.S76 I 

.063 and tha K-S probability la 6.1%. Flgura 9 glvaa tha 
usual (Y>P(Y)) probability curves for the power law alopa. Tha beat fit 
values for the paramacers are 


1st scan 


Y ■ 2.03 ± .18 

„ ,«+1.2. ..-7 .,-44 / .Y -1 „ -3 

K ■ (2 _) X 10 (10 ergs/s) Mpc 

" a O 


2nd scan 


Y - 2.07, 25 
K - (3.2 ± 2) 10 


(10^^ ergs/s)^ ^ 



The minimum luminosity Is now < 3.6 x 10^^ ergs/s (Abell 1060) In both first 
and second scan. The exponential fit Is again poorer# see Figure 10. 

2 . Active Galaxies 

1. Luminosity Function 

The Insert In Figure 11 represents the AML probability for the power 
law slope of the active galaxies differential luminosity function calculated 
from the 1st pass data. The best fit values for the power law parauneters are: 

Y « 2.75 ± *15 

K » (2.7 ± .15) X 10 ^ (10^^ ergs/s)^ ^Mpc ^ 

NGC 3227 Is the weakest source In the sample with 1.75 x 10^^ ergs/s and 
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XIIZW2 U th« brightest with i.3 x 10 ^^ ,rg»/». Figure 11 show* the binned 
representation (3 sourees/bin) » This result is similar to that of Boldt 
(1980) and Pye and Marwick (1979). The exponential form for the luminosity 
function is not acceptable as the probabilities are always less than 2%. 

The second pass sample is too small for a good determination of the 
luminosity function, however we find power law slopes steeper but consistent 
with the first pass ones 

li. h Lower Limit to the Active Galaxy Luminosity Function 

The active galaxies contribution to the cosmic X-ray background 
depends strongly on the lower luminosity limit of the luminosity function. 

The lower luminosity limit for which the function can represent the data, 
'^44HIN' ^ calculated by noting that the luminosity function must be 

consistent with the log N - log S observations. Namely, we can set a 
lower limit on by requiring that the number of active galaxies brighter 

than 1.25 R15 counts/sec expected from the luminosity function does not exceed 
the observed number plus 1 or 2 times the square root of the expected number. 

From eq (14.7.35) of Weinberg (1972), and assuming a power law 
luminosity function we have (for Y ^ 2.5 and Y ^ 3) 


NOS) » KA [L^'^”^] 


^MAX _-3/2 

Hiin 


B 

3-Y 



^MAX _-2, 


( 8 ) 


-9 -1 

where: S is the flux in ergs cm ^s 

K and Y are the parameters of the differential power law luminosity 
function in Mpe"^ (erg/sec)”^^ 

^MAX ^IIN lower limit of the luminosity function 

(actually NOS) depends strongly on and very weakly on L^;^) 
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ail th« luminoaitiaa ac« in args/i 
h •3.20 X lO*"^® 

B • 4.7 X 10"2® (assuming » 50 km/s/Npc) 

NOS) is tha total numbor of sources in the sky uncorrected 
for sky coverage. The second term of this equataion represents a first order 
cosmological correction to the Eucledian result. ^ ^ ^ 


( 1 ) Footnote I 

For L in units of lO^'* erg/sec equation (8) has constants 
A - 3.2 X lo"® 

B ■ 2.3 X 10"^ (Hq/ 50) (1+D where F is the spectral index of the source (here 
chosen to be .7) 


Assuming an average conversion factor of 2.17 x 10“^^ ergs cm"^s"^ per R15 
counts s“^ we find that the 1o lower limit on is 4 x 10^^ when y is 2.75 

and K is 2.68 x lO””^ (10^^ ergs/s)”^ Mpc""^ and varies between 5 and 15 x 

10^^ ergs/sec. 

In Table 4 we show the 1 and 20 limits on as a function of 

and y. We note that we have not included in Table 4 the possibility that all 
(or some) of the unidentified sources could be Seyfert galaxies. However, 
considering the distribution of identified sources with flux < 3 R15 cts/sec, 
we would expect, at most, 3 of these unidentified objects to be active 
galaxies. 

TABLE 4 

APPROXIMATE FOR VALUES OP AND y 

IWMt - 3 « 


- 1.5 » 10 
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Y 



Y 


2.6 

2.75 

2.9 

2.6 

2.75 

2.9 

1.5x10^2 

2.5x10^^ 

4.5x10^2 

1.5X10^2 

3.0x10^2 

4.5x10^2 

4.5x10^^ 

1.5x10^2 

2.5x10^2 

5.5x10^^ 

1.5x10^2 

3.0x10^2 


C • DiBcuBslon 

1. ClusterB of GalaxieB 

We note that our luminosity function for cluBtera of galaxiea 1 b 
very similar to the result of McKee et al# (1980). This indicates that, 
whatever selection effects are operating in making a X-ray or optically 
complete sample, they do not strongly bias the result. However there is a 
strong overlap in the individual objects between this sample and McKee's. The 
method we have used has allowed us in principle to discriminate between 
exponential and power law luminosity functions for clusters. It is somewhat 
surprising that a power law is favored, since it requires a change in form at 
low luminosities in order not to exceed the space density of all clusters 
(Bahcall 1979). However, the contribution of clusters to the diffuse X-ray 
background (DXRB) depends only weakly on the lower limit chosen. We do remind 
the reader that an exponential form is not excluded. Our data are not capable 
of rejecting the exponential form. They are also not capable of determining 
well the three constants in Bahcall* s (1979) suggested form of the luminosity 
function. 

Keeping in mind that the mean X-ray spectrum of clusters differs 
significantly from the diffuse X-ray background we shall, for historical 
reasons, compare the 2-10 keV volume emissivity of clusters to that of the 
diffuse X-ray background. For * 1/2, Hq = 50 km/sec/Mpc the 2-10 keV 
background has a volume emissivity of S 2.4 x 10^® erg/sec/Mpc^ . The 
i 

yi 


contribution of clustorn is 


j MIN w 1 X 10^® •r9«/8«c/Hpc^ 

%iAX 


{for max 10^® orgs/iec# ■ i x 10^^ •rg»/»«G, v>h«r« «• hav« usad 
the 1st pass cluster power law luminosity function without the Virgo 
cluster) . Therefore, in an average sense, clusters contribute ^ 4% of the 
2-10 kev background. (For a more accurate treatment of the problem which 
includes the effect of the spectral differences of clusters from the 
background see McKee et al. 19B0 and Marshall et al. 1980). We note that the 
present value agrees well with the estimate made by Marshall et al. (1980) of 
the maximum possible contribution of clusters if they were not to distort the 
thermal bromsstrahlung fit to the spectrum of the DXRB in tl>e 3-50 keV band. 
We note that the relatively soft spectra of clusters should result in an 
increase in their contribution to the DXRB in the Einstein Observatory energy 
range. 


2. Active Galaxies 

The luminosity function derived here is in reasonable agreement 
with those derived previously by Pye and Warwick (1979) and Boldt (1980) in 
both slope and normalization. Using a lower bound of 3.0 x 10^^ ergs/aec and 
a upper bound of 1.5 x 10“*® erg/sec for our luminosity function results in a 
volume emissivity of < 4.9 X 10^® ergs/sec Mpc® or a contribution of < 20% to 

the 2-10 keV DXRB. If the lower limit is Xo2 x 10^^ (see Table 4) the 
contribution to the DXRB is < 40%. In fact, in order not to exceed the DXRB 
the luminosity function of AGN's must flatten at L ^ 3 x 10^^ ergs/sec (De 
Zotti 1980). There is a strong Indication of such a flattening in the optical 
luminosity function (Huchra and Sargent 1973; Huchra 1977; Huchra 1980) at 
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My < •> 21.5 (Hq *50) •quivaiant to • optical bolomatrio lunlnoalty of ^ 1.2 x 
10^^ arga/aao. sinoa tha alopa of tha optical, lumlnoaity function* at hlghax 
luminoaitiaa* ia tha aamoi within arrora* (Huohra and Sargant 1973 > Waadman 
1979) aa tha X-ray function it ia tempting to aaaociata tha band in tha 
optical luminoaity function with tha band in tha X-ray function and tharafora 
dariva a ^ 35. Thia valua ia rathar larger than that found Ic^ examing 

individual objecta (Kriaa et al. 1980) Elvia et al 1978). Thia may be due to 
the fact that moat of the optical flux from low Iminoaity active galaxiaa 
doea not come from tha nucleua but from the atallar population. 

The total apace density Of X-ray emitting active galaxiaa in tha 
luminosity range 3 x 10^^ - 1.5 x 10^® la < 7 x 10"® Mpe"^ which ia < 1.5% of 
all galaxies of Mp < -19 (Huchra 1977). This compares to a space density of 
active galaxies of Mp < -19 of < 5 x 10“® Mpc“^ (Huchra 1977, 1980). It thus 
seems, to first order, thet all active galaxies of Mp < -19 emit X-rays at 
> 3 X 10^^ erga/aec. For a flat universe there ara (assuming no 
evolution) < 4 x 10^ X-ray emitting active galaxies with > 3 x 10^^ with 
3.5. 

We can also estimate, the number of sources per square degree expected 
in the Einstein deep survey if the luminosity function used in this paper does 
not evolve strongly in either slope or norm and that spectral effects, such as 
low energy absorption, are not important. With ■ .5, ■ 3 x 10^^ in 

the 2-10 keV band and, « 5 x 10”^^ ergs/cm^aec in the a^-IO kev band, 

(which corresponds to the Einstein "deep survey" limit for a a - 0.7 source we 
predict ^ 6 active galaxies per square degree and <1.3 clusters per square 
degree, compared to the 19 ± 8 total sources per square degree seen by the 
Einstein Observatory (Giacconi et al. 1979). DeZotti (1980) has performed a 
similar calculation and finds < 5 active galaxies per square degree for “ 
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9.1 X 10^^ and ■ 2.9 x 10^^ •xqu/»*e in tha 2-6 kaV band and aaauMiing 
that tha alopa of tha luminosity function is 2.5. Sinca most of tha objaots 
ara n«ar vm would axpact many of tha Binstain survay objaots to ba 
Sayfart galaxias of < 5 x 10^^ arg/sao and s ^ .20. This is a consaquanca 
of tha wall known fact that if tha luminosity function is staof>ar than 2.5, 
and t>arring strong avolution, whan ona looks at faintar objacts ona is looking 
primarily lowar in tha luminosity function rathar than at highar radshift 
objacts. 

h simple way to look at tha problem is to examine tha number of ob;jacts 
predicted our bast fit luminosity function which would have radshifts 
(z) ^ 0.5 and would have luminosities high enough to have bean included in the 
Einstein Deep Survey. (We shall use » .5 or 0 geometry for simplicity). 

For - 5 X 10"^^ ergs/cm^ sec in the 2-10 keV band and q^, * .5 that we 

predict < 1.4 X 10^ sources/ster due to active galaxies and < 1.4 xlo^ 
sources/ster due to clusters compared to tha 6.3±2.6 x 10^ sources/ster seen 
in the deep survey (Glacconi et al. 1979). We therefore predict that < 25 %-^^ 
of the sources in the deep survey are low (L 4 x 10^^) close by (z ^ .5) 
active galaxies or clusters of galaxies of luminosity > 1 x 10^^ erg/sec. 

That this was a likely situation was noted ly Fabian and Rees (1978). (If q^ 

- 0 the number of sources increases to ** 2.1 x 10^ sources/ster and the 
calculated contribution to the Einstein source counts to 35^^^%). 

Both the contribution of active galaxies to the DXRB and their 
contribution to the Einstein source counts depend sensitively on the lower 
limit, of the luminosity function used. It is possible that the 

luminosity where the flattening of the luminosity function takes place could 
be higher than our calcuated value if we allow a two slope model of the 
luminosity function rather than our simple single slope power law model with a 
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cutoff* Howavor our data art not good enough to conatrain auoh a modal* 
tharaforo atronqiy caution tha raaiar that thaaa raaulta ara modal daoandant 
and ahould ba treated an auoh* 

XX* CONCLUSIONS 

Wa have performed an all aky survey of x-ray aoiircaa complete to a 

limiting aanaitlvlty of 3.1 k 10“^^ arga/cm^ aac in tha 2-10 kaV band. Of tha 

85 detected sources only 7 remain without reasonable identifications. Tha 

log N- log S relation for axtragalactic sources is well fit by a Euclidean 

law ~ ■ 16.5 where S is in Rl5 ct/sac 

ds 

or ^ 2.2 X lo”^^ s”^’® (erg/cm^s)"^ sr"“* where s is in erg/cm^s in the 2-10 

keV band. This complete sample has allowed construction of luminosity 
functions based on a flux limited sample for clusters of galaxies and active 
galactic nuclei. These functions are well represented by power laws of slope 
2.05 and 2.75 respectively. The sample enables us to estimate that the 
luminosity function for active galaxies should flatten at L ^3 x 10^^ erg/sec 
in the 2-10 keV band. The cpace density of X-ray emitting active galaxies is 
approximately the same as that of optically selected Seyfert galaxies. ' 

Integration of the best fit luminosity functions indicates that clusters 
of galaxies contribute < 4% of the 2-10 keV diffuse X-ray background and 
active galactic nuclei < 20%. 'rne sum of these contributions is very similar 
to the 26±11% contribution due to resolved due to sources seen In the* Einstein 
deep survey. We also predict that many of the objects seen in the deep survey 
should be local# (z < 0.5), relatively low luminosity active galactic nuclei 
and clusters of galaxies. In order to determine more accurately the 
contribution of low luminosity active galaxies to the diffuse X-ray background 
one would have to sample the luminosity range over large solid 

angles. This would require a complete sky survey with ^ 30 times the 
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gansitivlty of the present one and a angular resolution ^ 20 times better. 
Such a survey would also extend the luminosity function up to luminosities 
of < 10^^ ergs/sec. We stress the importance of a conflate unbiased X-ray 
survey with good identifications in determining log N - log S and luminosity 
functions since there are various classes of sources of widely varying X-ray 
to optical luminosities. We feel that this strategy rather then deep 
observations over small solid angles will determine log N - log S and the 
luminosity functions most accurately for the local epoch since for a given 
observing time and fixed instrumental parameters the number of observed 
sources greater than some statistical limit is maximized when the solid angle 
is maximized at a given completeness level for a photon limited experiment. 
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riGURE ChPTIOMS 


Figu»« !• Th« conqpl«t«n«s8 leval of the pr« 80 nt tuev^y Vf «olt.ptio 
latituda. The diamond* are for the first pas* and the oro**e* for the second 
pa**> The lower histogram is the sky fraction in each ecliptic latitude bin 
(right hand scale). The centre of the diamonds and crosses is 5 time* the 
mean error for a source located tn that ecliptic latitude bin and the sire of 
the error bar is the standard deviation of this error. Since we truncate at 
1.25 R15 counts all of our sources at ecliptic latitude greater than 30^ lie 
well above the 50 level. We estimate that residual incompleteness of sources 
at levels less than 1.4 cts is less than 3 sources and zero source* greater 
than this limit. 

Figure 2. The error boxes for H0328-f025 and H09 17-074. The inner and outer 
boxes are the 90% confidence boxes as described in Marshall et al. 1979. The 
inner box assumes that the source was roughly constant during our period of 
observation • 


Figure 3. The distribution of the non-galactic sources detected in this 
survey in supergalactic coordinates. 

Figure 4. The probability distributions for K and a. The top panel shows 
the AML probabi^vity vs« a in the first pass data, the middle panel shows the 
AMTj probability vs. a in the second pass. The bottom panel shows the 68 and 
95% joint probability contour for K and o for the first pass data. The + 
marks the best fit. 


Figure $• The dltferentUl log N - log S diitrlbutlon Cor our ■•niple. The 
beet fit ie indiceted. The higheet flux point ie indiceted by a daehed oroea 
becauae ita upper flux bound ia not well defined, (let paaa data) 

Figure 6. The AML KolomogorovSmirnov teat diatribution for an a *> 2.5 
model. The SO and 95% probability .bounda are indicated, (lat paaa data) 

Figure 7. The ratio of the number of obaerved aourcea **obs to the number of 
expected aourcea Cor a >■ 2.5 log N - log S law. (lat paaa data) 

Figure 8a. The cluster of galaxies differential luminosity function for the 
first pass data. 

8b. The same information for the second pass data. The best fit 
power law models are indicated on both panels. 

Figure 9. the AML probability vs. Y the slope of the power law differential 
luminosity function for clusters of galaxies Cor the first and second passes 
including and excluding the Virgo cluster. 

Figure 10. Same as Figure 9 but for the exponential luminosity function.. 

Figure 11. The SeyCert galaxy luminosity function for the first pass data. 
The best fit power law differential model is indicated. The insert shows the 
AML probability vs. Y the slope of the luminosity function. 
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